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Abstract 

We propose the Kolmogorov stochasticity parameter, A for energy 
level spectra to classify quantum systems with corresponding classical 
dynamics ranging from integrable to chaotic. We also study the prob- 
ability distribution function (PDF) of A. Remarkably, the PDF of all 
the integrable systems studied here is the same and is found to be 
completely different from the PDF of chaotic systems. We also note 
that A„ for n energy levels scales as A^ ~ n~". Furthermore, with a, 
the stochastic probability (calculated from PDF) is seen to jump by 
about an order of magnitude as the systems turn chaotic. 
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1 Introduction 



Chaos in classical mechanics occurs as a result of loss of correlation between 
neighbouring initial conditions, leading to a very complex trajectory in phase 
space [1]. The faster this correlation is lost, more chaotic the system is 
and this is quantified by the Lyapunov exponent. However, in quantum 
mechanics, uncertainty principle does not allow a description in terms of 
trajectories which makes the notion of chaos in quantum systems to lose 
its meaning in the sense as it is used in classical mechanics. Moreover, 
unlike chaotic classical systems where one cannot obtain analytic expressions 
for classical trajectories, one can write (albeit rarely) such expressions for 
eigenfunctions [21 [31 SI [S] of some quantum systems which are classically 
chaotic. In this paper, we study the Kolmogorov stochasticity parameter for 
level sequences and several related aspects in attempt to capture order-chaos 
transition in quantum mechanics. 

We consider billiards which are dynamical systems consisting of a freely 
moving particle in an enclosure, reflecting specularly from the boundary in 
accordance with Snell law. These systems possess a wide range of dynamical 
behaviour depending upon the shape of the enclosure. For instance, circu- 
lar billiard enjoys two functionally independent constants of the motion, viz. 
angular momentum and energy, phase space trajectory resides on the surface 
of a two-dimensional torus. Such systems exhibit regular behaviour and are 
termed integrable systems if the dynamics is smooth everywhere [6]. How- 
ever, integrability is easily broken by breakdown of symmetry or by presence 
of singularities giving rise to diffraction effects [?]. A classical trajectory of 
an integrable system can be written in terms of Fourier modes built on pe- 
riods corresponding to the canonical angle variables, this is not possible for 
a chaotic system. Quantum studies [S] are concentrated on understanding 
the fluctuations in the level density in terms of invariant structures (e.g., 
periodic orbits) in phase space. However, for the billiards, the average level 
density is completely determined by the geometrical features like volume, 
area, perimeter, etc. owing to a well-known result of Weyl [S]. Thus, one 
concentrates on the fluctuations in the level density characterized by level 
correlations, which are eventually studied semiclassically in terms of periodic 
orbits by employing trace formulae. Perhaps the most popular measure is the 
nearest-neighbour level-spacing statistics [HI [10], P{S) giving the frequency 
of occurrence of certain spacing, 5* for a given spectrum of levels. For generic 
integrable systems, it was argued by Berry and Tabor that the spacing among 
nearest neighbours is distributed in a Poissonian manner [TT] . 

Statistical study of fluctuations of level density has a mathematical coun- 
terpart in statistics of eigenvalues of random matrices [I2] constructed within 
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certain well-defined symmetry constraints. If a physical system is invariant 
under rotation and time reversal invariance, then the spectral statistics could 
be compared to the corresponding random matrices that are real-symmetric. 
It is observed that the fluctuation properties of energy levels of quantized 
time-reversal invariant chaotic systems agree with that of orthogonal ensem- 
ble of random matrices exhibiting level repulsion as P{S) ~ S, for small 
S [ini [131 • However, a linear level repulsion is observed for quantized 
polygonal billiards [3 [HI [151 [IS] also even when these are non-integrable 
and non-chaotic [161 [13 [E] • This necessitates quest for a parameter which 
not only varies monotonically with underlying classical description but can 
also discriminate between chaoticity and non-integrability, independent of 
random matrix theory and semiclassical theories, however not in disagree- 
ment with them. 

One of the most significant advantages of this work is in its applicability to 
finite sequences of energy levels. In many physical situations, particularly in 
nuclear physics, the energy levels available are finite in number, of the order 
less than a hundred. In such situations, it becomes inappropriate to apply 
random matrix theory. We hope that this work will make a discussion of 
order-chaos transition in nuclei possible even for realistic systems with some 
tens of energy levels. Arnold |19| has employed the stochasticity parameter 
of Kolmogorov [20] to number sequences of length in tens successfully. This 
work is an adaptation of the work by Kolmogorov and Arnold to energy levels 
of quantum systems in particular, and to a sequence of modes supported by 
any physical system in general. 

The plan of the paper is as follows. In Section 2, we present the defini- 
tion of Kolmogorov stochasticity parameter, A. Based on semiclassical trace 
formulae, semiclassical expressions for the parameter are given for different 
systems and these are compared with the numerical calculations. Putting all 
these expressions together, we show that A scales with the length of the set 
of energy levels. In Section 3, we recapitulate Kolmogorov statistic briefly 
and present some of the applications of this for number sequences, presented 
recently by Arnold [19]. Subsequently, we present probability distribution 
functions (PDF) of A for energy levels of different dynamical systems. We 
also present results for zeros of Riemann zeta functions, and, for Gaussian 
unitary random matrix ensemble. Finally, we conclude by showing that with 
these PDFs, we can visualize an order-chaos transition in spectral properties 
of dynamical systems. 
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2 Kolmogorov stochasticity parameter 



2.1 Definition 

Consider a quantum system Si which possesses energy levels corresponding 
to its bound states, Ei < E2 < ■ ■ ■ En. Cumulative density of levels is defined 



as 



Nn{E) = Y,9^Q{E-E,) (1) 

i=l 

where Qi denotes the degree of degeneracy of the level, E^. The average 
cumulative density of levels can be obtained by various methods for quantum 
systems. To be specific, we take simple systems that are illustrative and very 
popular in studies in quantum chaos - the two-dimensional quantum billiards 
and quantized maps. Quantum billiard is a system consisting of a point 
particle of mass m enclosed in a box of a certain shape, the solutions of 
the Schrodinger equation are assumed to exist and the eigenvalues {E^} are 
obtained. For such systems, the average cumulative density of levels, Nq{E) 
is given by Weyl formula [2T]. 

In analogy with the definition given by Kolmogorov [22j, let us define 
stochasticity parameter, 

iN^jE) - NoiE)\ 

A„ = sup 1= . (2) 

E -Jn 



In the limit of n — )■ 00, the difference Nn{E) — No{E) will be just the oscil- 
lating part of the density of levels for which we can employ trace formulae 
[S] whenever possible. For the eigenvalues of a Laplacian (sometimes related 
to number theory), there are trace formulae which can be employed one 
of the well-known being the Selberg trace formula. This further helps us to 
show that even finite sequence length of size 5000 reproduces the behaviour 
obtained from trace formulae. In the Appendix, we have collected the nu- 
merical calculations of A„ and compared them with the results based on trace 
formulae. As can be seen there, these results agree well with each other for 
various billiard systems and maps. 

We would like to recall that statistical investigations of Nn{E) — Nq{E) 
have been carried out in the past, albeit with a different normalization, as a 
measure of quantum chaos in spectra [23j. In these studies, the measure was 
taken as 

Wi.) ^ (3) 



4 



where Aoo{E) is the hmit of the spectral rigidity, A-^{L; E) as L — )■ oo. It 
was shown that W{x) possesses a hmit distribution with zero mean and unit 
variance for bound, conservative and scahng systems. For classicaUy (inte- 
grable) chaotic systems, the distribution was expected to be (non-) Gaussian. 
These predictions were verified in a statistical analysis of eigenmodes of 
superconducting microwave billiards [23]. For a Bunimovich stadium, the 
three-dimensional Sinai billiard, and for Limacon billiards, a Gaussian dis- 
tribution was found. In contrast, for a two-dimensional circular billiard, a 
non-Gaussian distribution was found. 

Practically, to calculate A„ the numerical procedure followed is sketched 
as follows: (i) sort the given sequence; (ii) find out how many eigenvalues 
are less than or equal to E, from this subtract the expected number of such 
eigenvalues (calculated from Weyl's formula) and divide by square root of 
this number; (iii) vary E and then take maximum value of number calculated 
from the last step, this is A„ for fixed sequence of length n. 

2.2 Semiclassical expressions 

Gutzwiller discovered the connection between density of energy levels of a 
quantum system and the periodic orbits of the underlying classical system 
[25] . The ensuing trace formulae have been found for quantum systems with 
underlying classical dynamics from regular to fully chaotic. For the harmonic 
oscillator potentials, the trace formula turns out to be analytically exact [26]. 
In the following, we use the well-known trace formulae to get the Kolmogorov 
stochasticity parameter, X{E) and compare the expressions with numerical 
results. In the subsequent Sections, we will employ these to study probability 
density distribution functions of A(-E'). Thus, the agreement between the nu- 
merical and semiclassical results allow us to decide the length of the sequence 
required in studies in latter Sections. As seen in Fig. [H the agreement of 
numerical and semiclassical results is very good. 

2.3 Scaling 

Studying A„ for various systems discussed above as a function of In n (Figj2|) 
reveals an exponential fall for Riemann zeros and weaker exponential curves 
for other systems. Indeed, on comparing each of them with exp(— alnn), 
we extract a monotonically changing index, a varying from 0.15 for circular 
billiard to 0.45 for Riemann zeros with a value 0.34 for the chaotic standard 
map. The values of this index are listed in Table I. These show an increasing 
trend as the system becomes classically chaotic. 
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(a) rectangular billiard 
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(b) circular billiard of unit radius 
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(c) triangular billiard of side length 
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(d) zeros of Riemann zeta func- 
tion, the values of a is taken as \ 



Figure 1: Stochasticity parameter calculated using the exact expression and 
trace formula, the agreement is good. 




Figure 2: Stochasticity parameter for all the systems considered are displayed 
here as a function of logarithm of the length of the sequence, In n. This brings 
out an interesting observation, A scales as n~". This leads us to finding the 
best fitted values for the index, a. As the systems become increasingly 
stochastic, the value of the index increases. 
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We may understand the scaling exponent a heuristically in the following 
manner. For integrable systems, like rectangle the trace formula contains 
Bessel function and a trigonometric function for i?-dependence apart from 
E~'^/'^ dependence in the denominator. For large Ji{z) ~ ^y2/{'Kz) cos [2 — 
3'7r/4]. After a little simplification one will see that 

A„ ~ sup \[E-'/\^{^/E) + MVe) + MVe)] 

E 

- E-'%,{Ve) + g2{VE)]\. (4) 

where /j's and gi^s are trigonometric functions (sin, cos) of \/E. This leads to 
an scaling with a maximal exponent of about 0.25 in asymptotic limit with 
an oscillatory dependence with an argument containing \/E. 

For the Riemann Zeta function it is just the denominator modified by the 
oscillatory behavior of oscillatory function which, to a good approximation 
may be written as t^'^. 

Note that the above heuristic argument give a ratio of A for the zeta 
function (=0.6) to that for rectangle (=0.25) is 2.4 while the numerically 
found ratio 0.4451/0.18 is 2.47 - they are close. 

From [23], one may expect scaling behaviour because the systems consid- 
ered here are scaling systems in the sense discussed above. 



3 Probability distribution functions 
3.1 Kolmogorov statistic 

Since the average level density used in the definition of A is a continuous 
function, it renders a continuum character to A itself which in turn makes it 
possible to obtain the cumulative probability distribution function, $ (A) of 
variables, A. $(A) gives the probability of the stochasticity parameter to have 
a value A < A. In the case of independent identically distributed random 
variables, this distribution is indeed a limiting distribution of Kolmogorov 



$(A) = ^ (-l)'=exp(-2A;2A2) (5) 
fc=— 00 

It grows from $(0) = to $(+00) = 1. The mean value of A is \/tx jl log 2 ~ 
0.87. This distribution has been used for comparing the degree of stochastic- 
ity in geometric and arithmetic progressions, among other examples E]. 
Larger stochastic probability is obtained for a geometric progression than an 
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arithmetic one, with only fifteen numbers. The stochastic probability of a 
sequence is quantified by the value $(A'') where A° is obtained according to 
the definition above. 

3.2 Extensions of Kolmogorov distribution for dynam- 
ical systems 

As suggested by Arnold [31] , it will be very interesting to extend Kolmogorov 
formula to the more general situation for random variable having continuous 
cumulative distribution function. For the dynamical systems considered here, 
we have such a continuous distribution function. This gives us an opportunity 
to propose some of the extensions of Kolmogorov formula. We present results 
from our numerical investigations as well as analytic calculations based on 
trace formulae. 

The numerical evaluation of a probability distribution function is sum- 
marized in the following steps: 

• for an ordered sequence of length n, we can have n-subsequences of 
length 1, 2, ... n, always starting from the lowest element and increasing 
by one element; 

• for each subsequence we calculate A; 

• this way, we end-up with n. As and a histogram of these As gives the 
probability distribution function (PDF) of A for the sequence. 

3.2.1 Integrable billiards 

For the integrable billiards like circle, rectangle, and equilateral triangle, 
the values of the stochasticity parameter, A are respectively 0.19, 0.13, and 
0.12, all comparable; these values are stable with the length of the sequence. 
Probability distribution function of A for the integrable systems are very close 
to each other (Fig. [3]) fitting the function, ^integratiei-^) ~ Ciexp(— A'') (for 
details, see the Table. The value of $ at A gives the stochastic probability 
of the level sequence, in line with the measure of stochasticity for number 
sequences developed by Arnold [19]. 

3.2.2 Non-integrable systems 

For chaotic systems, as mentioned above, spectral fiuctuations are found to 
be in close agreement with random matrix theory. A convincing argument 
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Circle billiard 




Rectangle billiard 




Triangle billiard 




\\ ' ' 


^v' - V 



Stochasticity parameter (X^) 

Figure 3: The probability distribution of the stochasticity parameter is shown 
for three integrable bilhards, viz. circle (red), rectangle (blue), and equilat- 
eral triangle (black). For these billiards, the distribution can be fitted to the 
same functional form, Ci exp(— A^)), with (ci, f3) taking values (4.344, 2.927), 
(4.163, 3.196), and (4.703, 3.826) for circle, rectangle, and equilateral trian- 
gle respectively. Value of the stochastic probability at A measures quantum 
stochasticity. 
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Figure 4: The probability distribution of the stochasticity parameter is shown 
for standard map, GUE, and Riemann zeros. For these systems, the distri- 
bution can be fitted to the same functional form, (c+ A)~^, with (c, 7) taking 
values (0.1688, 3.632), (0.1932, 4.097), and (0.2103, 4.397) for standard map, 
GUE, and Riemann zeros respectively. Value of the stochastic probability at 
A measures quantum stochasticity. 
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for this statement has been given recently[3n] on the basis of Gutzwiller peri- 
odic orbit theory [Bl[in]- For the quantized standard map, Xstandard for 5000 
levels is 0.0486. The function that fits the probability distribution of A for 
standard map, Riemann zeros, and eigenvalues of Gaussian Unitary Ensem- 
ble are all close to each other, represented by [c + \)~'^ where 7 is about 3.6 - 
4.4 (see Tabl^. We find the measure of stochasticity, ^{Xstandard) = 0.4862. 
Further, on considering zeros of Riemann zeta function, = 0.0315 and 
$(A(^) = 0.3775. These are more than twice compared to the values for in- 
tegrable case. Since it is very well-known that the fluctuation properties of 
Riemann zeros agree with those of the eigenvalues of unitary ensemble of 
random matrices, we have also calculated the stochasticity parameter and 
its distribution for them. We find Xgue = 0.035, fitting function is identical 
and ^{Xgue) = 0.4028. For the 7r/3-rhombus billiard, classical phase space 
surface is topologically equivalent to a sphere with two handles [IZIE], Lya- 
punov exponent is zero. The energy levels of this almost integrable system 
entails an intermediate value for the stochasticity parameter. The probability 
distribution function also shows a concurrence with the family of integrable 
billiards (see Table [3]) . 

3.3 Order-chaos transition 

From the nature of probability distribution functions, we observe that the 
derivative of $ at A calculated for different systems jumps by six to eight 
times (Table [1]), thus behaving like an order parameter for the order-chaos 
transition. The distinction among the cumulative distribution functions of 
integrable, chaotic and Kolmogorov distribution is apparent in Fig. [51 These 
results point at certain interesting conclusions: (i) stochasticity parameter 
can characterize the level sequences coming from random matrices, quantum 
systems, number theory etc., (ii) probability density distributions are dis- 
tinctly different for integrable and chaotic quantum systems, (iii) we observe 
the confluence of the statistical behaviour of Riemann zeros and unitary 
ensemble in terms of A where we have not imposed any symmetry-related 
constraints as in random matrix theory. In a nutshell, we show (FiglH]) the 
value of probability distribution function $' for all systems with a. 

3.4 Connections with number theory 

Finally, we would like to present the precise connection of our ideas with a 
well-known number-theoretic problem of counting the lattice points in a circle 
of radius, R, first studied by Gauss^38j. The difference of this number with 
the area of the circle, known as the error term with proper normalization. 
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Figure 5: A comparison among CDF for three different class of dynamical 
systems 
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Figure 6: Two measures of stochasticity are being plotted here against each 
other. For calculating the PDF at stochasticity parameter for a specific 
system, we have normalized the distributions. As the index a increases, we 
see that the probability distribution of the stochasticity parameter undergoes 
a jump by an order of magnitude. This — a graph is a defining feature of 
quantum chaos. 



11 



Table 1: From the variation of stochasticity parameter, A with the length 
of sequence, n (Fig. [2]), we extracted a scaling exponent a. Employing 
the normalized cumulative distribution functions (CDF), we have also noted 
that the value of these are considerably larger for chaotic systems than for the 
integrable cases. As we go down the columns, we obtain larger values for the 
exponent and values of $' increase too. However, it is the PDF (rightmost 
column) that jumps up to eight times which is enormous for the fact that all 
the distributions are normalized. 



System 


a 


$ 


$' 


Circle 


0.1509 


0.2126 


1.11239 


Rectangle 


0.1757 


0.1523 


1.11465 


Triangle 


0.1808 


0.1340 


1.10569 


Rhombus 


0.2519 


0.0479 


1.09771 


Standard 


0.3409 


0.4862 


6.22019 


Map 








CUE 


0.3874 


0.4028 


8.10375 


Riemann 


0.4451 


0.3775 


8.74426 


zeros 









has been investigated extensively [39]. A connection of this problem with 
quantized integrable systems in two dimensions which have two quantum 
numbers has been studied in the past [10], however a measure of stochasticity 
was not discussed. The emergence of probability distribution functions of 
A for integrable and chaotic systems has been shown here. For cases that 
cannot be cast in terms of a problem with lattice points, our discussion opens 
up an enormous mathematical challenge to develop arguments which would 
probably result from a generalization of the classic work of Heath-Brown . 

4 Concluding remarks 

As the nature of classical behaviour of a system becomes more complicated, 
we have discovered the existence of a scaling exponent a of stochasticity 
parameter with the number of levels which follows the classical trend. Our 
studies are based on a number of systems displaying a variety of dynamical 
behaviour which have been shown to belong to distinct classes, not only in 
terms of this exponent but also in terms of the probability density distribution 
functions, $' of the stochasticity parameter. The values of $' show a jump 
by a factor of six to eight as the system becomes chaotic. 
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Table 2: Comparison of Stochasticity parameter for three different class of 
dynamical systems 



General model:/(x) = o * ex.p{—x^) + 0.0025 


Rectangle 


Circle 


Triangle 


Coefficients (with 91 
% confidence bounds) 
a=4.163 (3.575,4.75; 
b=3.196 (2.202,4.19) 
Goodness of fit 
SSE: 2.33 R-square 
0.8721 A5000 = 0.1365 
$'(A) = 1.1146 




Coefficients (with 95% 
confidence bounds) : 
a=4.344 (3.907, 
4.781) b=2.927 (2.304, 
3.551) Goodness of fit: 
SSE: 1.509 R-square: 
0.9227 A5000 = 0.1900; 
$'(A) = 1.1124 


Coefficients (with 95% 
confidence bounds): 
a=4.703 (4.101, 
5.305) b=3.826 (2.719, 
4.933) Goodness of fit: 
SSE: 1.84 R-square: 
0.9023 A5000 = 0.1212; 
$'(A) = 1.1057 


General nio(lel:/'(,7-) = (« + ,/•) 


Standard Map 


CUE 


Riemann zeros 


Coefficients (with 95% 
confidence bounds) 
a = 0.1688 (0.1606 
0.1771) b = 3.63S 
(3.463, 3.8) 
Goodness of fit: SSE 
0.01436 R-squarc 
0.9972 A5000 = 0.0486 
$'(A) = 6.2202 




Coefficients (with 
95 % confidence 
bounds): a = 0.1932 
(0.1907, 0.1958) b = 
4.097 (4.02, 4.173) 
Goodness of fit: SSE: 
0.000886 R-squarc: 
0.9997 A5000 = 0.035; 
$'(A) = 8.1038 


Coefficients (with 
95 % confidence 
bounds): a = 0.2103 
(0.2091, 0.2115) b = 
4.397 (4.336, 4.459) 
Goodness of fit: SSE: 
0.0001541 R-squarc: 
0.9999 A5000 = 0.0315; 
$'(A) = 8.7443 



Table 3: Fitting parameter for Rhombus 



General model: 

f(x) = a* exp{-x.^) + 0.0025 
Coefficients (with 95 % confidence bounds): 
a = 3.626 (3.534, 3.718) 
b = 4.367 (4.092, 4.642) 
Goodness of fit: 

SSE: 0.02527; R-square: 0.9972; 
A5000 = 0.0436; $XA) = 1.0977 
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In his Kolmogorov lecture Sinai explained that Renormalization 

Group method is connected to scaling ideas and limit theorems in probabil- 
ity theory. With the scaling property of A„ and existence of distinct PDFs, 
we may speculate a renormalization group argument to formally explain the 
transition noted in this work. 
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Appendix 

Here we collect the trace formulae [H] for different classes of dynamical 
systems and present the calculations of the stochsticity parameter based on 
them. The point is to make a comparison with the numerical results and 
we shall see that with 5000 energy levels, we have a good agreement (Fig. 
1). We have chosen the following examples to guide us to the length of the 
sequence of levels required to arrive at reliable conclusions: (i) rectangle bil- 
liard - a separable, integrable system where the quantization condition is an 
algebraic relation; (ii) circular billiard - a separable, integrable system where 
the quantization condition is a transcendental equation; (iii) equilateral tri- 
angle billiard - a non-separable, integrable system where the quantization 
condition is an algebraic relation; (iv) Riemann zeta function - nontrivial 
Riemann zeros are being considered here, the interest being due to their pos- 
sible connection with unitary ensemble of random matrices which, in turn, 
share spectral fluctuations with some well-known quantum systems which 
are classically chaotic. 
Rectangle billiard 

Let us consider a rectangle billiard where the sides of the rectangle are 
ai and 02 for which the density of energy levels (derivative of the cumulative 
density) is given by the trace formula [8]: 



'"Vir) 



i=l 

(00,00) 

maia2 

(Mi,M2)=(-oo,-oo) 

a,' / Im 



E£^Vf^ E cos (2Ma.v^M) . (6) 



1=1,2 Af=-oo 



where Jq{x) is the cylindrical Bessel function [31] and action, SmiMz is given 
by 



SmiM2 = V2mELM,M2, Lm,M2 = \ Mfaj + M^aj. (7) 
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With this, putting ^ = 1 (for all the systems 

("l."2) 

4vr ..-f^. , Lmi,M2 



N^{E)-No{E) = ^ T^M^Lmm 

(Ml, Ma) 



+7? E ^^.l^i«..ol 

Mi=l ^"1'*^ 

aia2 -A ^ 
-— ^ — (^sin(2Mai/E) + sin(2Ma2/E))(8) 

M=l 

where Ji{x) is cylindrical Bessel function of order one. 

For the denominator of ([2]), we will use Weyl's staircase function for the 
energy levels of rectangle. With the corner corrections included, this is read 

as 

n{E<E,) = — + -. (9) 

where area and perimeter is taken vr, 2(1 + vr) respectively. 

Using first 5000 energy levels of rectangle, we calculated Kolmogorov's 
stochasticity parameter ([2]) as a function of length of energy level sequence 
and compared the same when is calculated using and (jH]). The agree- 
ment is shown in FiglTal 
Circular billiard 

For a circular billiard of radius R, the treatment in terms of periodic 
orbits can be seen in the work of Balian and Bloch [29]. However, the trace 
formula was given by Bogachek and Gogadze [3D]. The oscillatory part of 
the density of eigenvalues is 



gosc{E) = ^.l — ^y^y^ fvw ^^sin$„^, (10) 

" 10=1 v=2'w 



where 



fx TXUJ 

t; 1 'Pvw ; -^tito 2'1'i? sin 0^^, 



P 



V2mE, ^^y, = — 3t;- + — , 

n 2 4 

fyw = 1 ioT V = 2w 

= 2 for v> 2w. (11) 
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The average cumulative density of energy, Eq for a circle of area a, cir- 
cumference, 7 is given by [2T] 



n{E < En) = — . 

^ - 47r 47r 6 



:i2) 



Nr,{E)-No{E) 



^= > / dE' — 



4:wRV2mE' - Sttw + ^ 



{2mE')-i 



AmR^ \ ^ \ ^ 1 3 /7rw\ 
, — > > ^sin2 — 

" ui=lii>2ui * 



B sin |2t;i?sin^V2mE' - 3t;f + ^ 
dE'- 





The integral involved in this expression is 



{2mE')\ 



■ (13) 



I{E) = / sin(7A/x + 5). (14) 

Let ^/x = {y — so dx/{2y/x) = dy/'j. Then, 1(E) becomes 



HE) 



tVE+S 



dyz^smy{y - 6) 



1/2 



7 



sm{z + 6)z^^'^ [z = y — 5) 



7 



3/2 



cos(7V-c/ + o) + A/ ^ coso IresnelC 



7 



'27V^ 



TT 



-FresnelS 



/27/E 




V . 


sin 5 J 



(15) 



In f[T5|) . the terms containing FresnelC and FresnelS terms contributes at the 
third decimal place which for all the numerical comparison purposes we can 
ignore. The agreement is good and given in FigJlbl 
Equilateral triangular billiard 
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Equilateral triangular billiard represents a class of dynamical systems 
which are integrable but not separable. For equilateral triangular billiard, 
the eigenvalue spectrum is given by 

E(rn, n) = {jn? + 11? — mnj 

16 h'^Tv^ 

^^-^2^? " = 1.2.3.... (16) 
n = 1,2, 3, ...(m > 2n). 

All the eigenvalues are doubly degenerate except (m = 2n) which is singly 
degenerate. Taking these in consideration the level density is given by 

3V3Ea M„tt=-oo 2^EEa ^f^_^ 



where Lm^M2 = y^i^f + M| + MiM2)L. 
The analogue of ([HD and for equilateral triangular billiards are 



Nr,{E)-No{E) 



TT 



f: ^MVELo,M2] 



TT 



3\/3Ea Mfji2=i 



TT 



00 —1 



7 Ji[VELm,,M21 



2 °° 1 

sin(v/ELM) 



"^En I En 1 , , 

In our calculation of stochasticity parameter for this system Ea = tt^. 
A comparison similar to rectangular billiard between the stochasticity pa- 
rameter calculated exactly using 5000 energy levels of equilateral triangular 
billiard and using trace formula is presented in FigJTcl 

Zeros of Riemann zeta function 



18 



Riemann zeta function is a very important special function [31J with 
many interesting connections with number theory. It is defined in terms 
of a Dirichlet series or a product form as follows: 



E- 

n=l 

n 

primes, p 



1 - 



p, 



-1 



(20) 



According to the Riemann hypothesis, all the nontrivial zeros lie on the 
critical line, given by s„ = ^ + itn- The product formula is convergent only 
for 0" > 1, but we employ it to write an analogue of Gutzwiller trace formula 
[32]. Taking logarithm of the product form above and using log(l — x) = 



exp(— ifct log p) 

f^pkcr 



(21) 



Taking imaginary part of both sides, and differentiating with respect to t, we 
obtain for fixed a > 1: 



^o- = -^EE^^°^(^^l°gP)• 



p k ^ 



The oscillating part of the cumulative density will be 



TT — ' — ' fcn^ 

p k 



(22) 



(23) 



For the denominator of ([2]) 
Int 



n 



t , t t 
— log 

271 ^ 27r 27r 



(24) 



Finally, for the zeros of Riemann zeta functions on the critical line, 



sup 

t 



IN'' I 

I osc I 



sup \ - 

t V vr 



oo 



sin [(A; logp)t] 



(25) 



This expression should be doomed as its convergence is in the domain where 
there are no nontrivial zeros. However, as noted in [33], trace formula con- 
tains information about the zeros on the critical line, with the shortest orbits 
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(corresponding to the smallest primes) contributing the most P?]. We have 
compared the values of stochasticity parameter using numerical values of the 
zeros (www.dtc.umn.edu/~odlyzko/zeta_tables/) and the expression (125|) in 
FigJldl For calculation we have utilized first 5000 only. This is done to 
compare the results from various systems for same sequence length. 
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